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We consider a gas of bosons in a trichromatic optical lattice at finite temperatures. As the 
amplitude of the secondary lattice grows, the single-particles eigenstates become localized. We 
calculate the canonical partition function using exact methods for the noninteracting and strongly 
interacting limit and analyze the statistical properties of the superfluid phase, localized phase and 
the strongly interacting gas. We show that those phases may be distinguished in experiment using 
off-resonant light scattering. 
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I. INTRODUCTION 

In 1958 Anderson proposed a mechanism that ex- 
plained the absence of conductance in certain types of 
media |1 . By considering the propagation of mat- 
ter waves through a medium containing randomly dis- 
tributed impurities, he discovered that under some con- 
ditions the single-particle states become exponentially lo- 
calized and the propagation is blocked. Anderson argued 
that the effect occurs due to destructive interference of 
waves, therefore it can be observed both for quantum 
matter waves and for classical waves. So far Anderson 
localization (AL) has been reported for acoustic waves 
[2], microwaves [3J, and light [3]. It has been also real- 
ized in ultracold atomic gases confined in optical random 
potentials [5J [B]. In comparison to other physical sys- 
tems, ultracold atoms offer unique possibilities of control 
of both the shape of the external potential and atomic 
interactions [7J. Disordered potentials can be generated 
either by use of the laser-speckle pattern [SJ [3] or a two- 
color optical lattice PJJ [TJJ . In the former case the po- 
tential is truly random, while in the latter the system is 
quasi-periodic. 

Localization properties depend strongly on the form of 
the potential and the dimensionality of the system. It 
is known that in one dimension AL is present even for 
arbitrarily small random potential [T2TI14) . However, for 
a pseudo-random potential, such as the two-color optical 
lattice, there is a transition to localized states only at 
certain disorder strength |15j . 

One of the most important challenges is to under- 
stand the interplay between the disorder and interac- 
tions. Adding the interactions may lead to novel quan- 
tum phases, such as the Bose glass phase, which ap- 
pears in addition to superfluid and Mott insulator phases 
[To] HT] , The phase diagram of such systems can be 
complicated and hard to obtain theoretically (see e.g. 
[T5lll9| ). Another challenge is to include the role of finite 
temperature, which excites the particles and suppresses 
localization. 

Detecting and studying the correlation properties of 
the quantum phases in experiment is not easy as well. 
The usual detection scheme is to switch the trap off 
and let the atoms expand for a certain time, measuring 



the interference pattern 20J . However, the information 
provided by this method is very limited [21 , in partic- 
ular in case of disordered lattices. Therefore new ex- 
perimental schemes has been developed, e.g. based on 
the noise correlation measurements [22] [23] , or observing 
the far-field pattern in coherent light-scattering |24j . So 
far methods based on atom-light interactions have been 
proposed for several various situations, such as detection 
of the Bose condensed phase [25], condensate fluctua- 
tions [25] , superfluidity in Fermi gases [27] and for detec- 
tion of quantum phases in optical lattices [28H3T] . The 
spectrum of weak and far-detuned light carries informa- 
tion about the static structure factor of the investigated 
system [21] [29] [31], which contains information about 
the density-density correlations and the energy spectrum. 
The structure factor is also accessible via Bragg diffrac- 
tion [2TJ [32] • The phases of strongly correlated systems 
can be also detected using quantum-noise-limited polar- 
ization spectroscopy [33] . 

In this work we investigate the finite temperature prop- 
erties of a one-dimensional Bose gas in a bichromatic op- 
tical lattice. We focus on two limits when the partition 
function can be found exactly. One is the noninteract- 
ing gas described by the Aubry- Andre Hamiltonian and 
the other is strongly interacting gas with negligible inter- 
well tunneling. We calculate the mean, fluctuations and 
correlations of occupation numbers in the lattice wells. 
In addition, for the ideal gas we examine the condensate 
fraction and its fluctuations. In the second part of the 
paper we analyze the angular distribution of light scat- 
tered from the ideal and strongly interacting gas. We 
show that similarly to regular optical lattice the light 
can be used to discriminate different phases existing in 
disordered systems. 

Our paper is organized as follows. In section II we 
introduce the Aubry-Andre model, review its single- 
particle properties and analyze the energy spectrum. Sec- 
tion III is devoted to the statistical properties of the gas 
in two-color optical lattice and the impact of temper- 
ature on the localization properties. We consider the 
strongly interacting gas in section IV. Section V studies 
the properties of the off-resonant light scattered on non- 
interacting and strongly interacting gases in disordered 
potentials systems and discusses distinguishability of dif- 
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ferent phases with this method. Section VI presents the 
conclusions and three appendices give some technical de- 
tails on exact calculations of the statistical quantities in 
the canonical ensemble. 



II. BOSE-HUBBARD MODEL FOR 
DISORDERED SYSTEM 

We consider a ID Bose gas in the optical lattice in the 
presence of a weak disorder. We assume that the gas 
is sufficiently cold and its dynamics takes place only in 
the lowest Bloch band. In such a case the dynamics is 
governed by the Bose-Hubbard Hamiltonian [34] 

H = - J Y 9i9j + j Y n ^ Ul ~ ^ + Y Ul£l w 
(id) 1 i 

with an additional term describing the on-site energies 
due to the disorder [TTJ [TB] . Here, <?, and g\ are the an- 
nihilation and creation operators of bosons at lattice site 
i, respectively, ni — g\gi is the particle number operator 
at site i and J and U are the energy scales corresponding 
to the tunneling between wells and the on-site interac- 
tion. In our approach we consider disorder induced by a 
bichromatic optical lattice potential [15] 

V(x) = s\E r i sm 2 {kix) + S2-EV2 sin 2 (fc2^ + </>). (2) 

Here, /cj = 2ir/\i are the wave numbers of the light beams 
creating the standing wave, Ei = h 2 /2mA 2 are the recoil 
energies, Sj are the heights of the two lattice potentials in 
units of recoil energies, cj> is a relative phase between two 
laser beams and m is the atom mass. We will denote the 
ratio of the wave numbers ki/k% as /3. For s\E r i 3> $iE r i 
the first lattice generates the periodic structure of the 
potential, while the second lattice generates weak quasi- 
periodic modulation of the potential wells. In this case 

= AcOs(27T/3i + 2(f)). (3) 

Here, A is the measure of the disorder strength. It 
can be expressed in terms of the Wannier states w(x) 
localized in the wells of the first lattice [35] : A = 

£ rl -^/decos(2/90MO| a - 

The case of an ideal gas U = at zero temperature 
has been extensively studied in the literature 1151 1351 136] . 
The disorder introduced in this Hamiltonian is pseudo- 
random, and it has been shown that even in one dimen- 
sion the eigenstates of the single-particle Hamiltonian are 
not localized for low A, in contrast to the standard An- 
derson localization [14] . Instead, if (3 is an irrational 
Diophantine number, there is a transition from extended 
to localized states. In particular case of /3 = (V5 — l)/2 
the transition occurs at A — 2 J, which is a self dual 
point [151 137] . In practice, the system of atoms in an op- 
tical lattice has a finite size, therefore it is sufficient that 
ft is a rational number and the periodicity of the on-site 
energy modulation e n is larger than the system size. 



III. IDEAL GAS 

In this section we consider the statistical properties of 
an ideal gas confined in a quasi-periodic potential. We 
start by studying the single particle properties, then we 
study the statistics of the gas at finite temperatures. 

A. Single-particle states 

In case of an ideal gas (U — 0) the elementary excita- 
tions of the Hamiltonian ([!]) 

H = -J Yah 1 + A^n fc cos(27r/3fc) (4) 

(k,l) k 

have been already studied by Aubry and Andre [TS]. As 
the Hamiltonian is quadratic it can be easily diagonalized 
in the basis of states describing atoms localized in a sin- 
gle potential well. In this way the single-particle states 
\^n) and corresponding single-particle energies e„ can be 
expressed as 

F|* n ) = e n |*„), (5) 

where |0) denote the vacuum state, and c" are expansion 
coefficients. We calculate the energy spectrum and eigen- 
states of the model numerically. In original formulation 
of Aubry- Andre model (3 is irrational and the system is 
quasiperiodic. In this case it is impossible to use periodic 
boundary conditions. In practice, however, it is sufficient 
to use a rational j3 with sufficiently large system so that 
the periodicity of the energy modulation e m is equal to 
the system size. A practical way to do that is to use 
Fibonacci numbers fi, setting f3 — f n -i/fn and M (the 
number of lattice sites) to f n . Most of the numerical 
results in this work are obtained using M = 144 and 
f3 = 89/144, which is a sufficiently good approximation. 

Figure [l] shows the energy spectrum for different 
strengths of the disorder. When A = the energy 
spectrum is just a lowest Bloch band, like in the Bose- 
Hubbard model. When the disorder increases, some en- 
ergy levels tend to form groups, separated from the others 
by energy gaps. The effect is strongest for A = 2 J, for 
even larger A the spectrum becomes more regular again. 
The behavior of the ground and the first excited state 
for different values of A is shown on Figs. [2] and [3] The 
ground state for A = is almost uniformly distributed 
over the whole lattice. When A gets larger the proba- 
bility distribution becomes nonuniform and some lattice 
sites are favored. Finally for A = 2 the ground state be- 
come exponentially localized in a single site. In contrast, 
for the same value of A the first excited state exhibits 
two maxima localized in two distant lattice wells. For 
even larger A, however, the first excited state becomes 
localized in a single lattice well. We have verified that a 
similar behavior occurs for higher excited states. 
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FIG. 1: Energy spectrum of the Aubry- Andre Hamiltonian 
Q for M — 610 lattice sites and for different values of the 
disorder amplitude A. 
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FIG. 2: Probability distribution |cf| 2 of finding atoms in dif- 
ferent lattice sites for atoms in the ground state and for vari- 
ous disorder strengths A. We use logarithmic scale for y-axis 
to show exponential localization. 



B. Statistical properties at finite temperatures 

We now examine the properties of noninteracting 
gas in a Dichromatic lattice at finite temperature. As 
the grand-canonical ensemble predicts unphysically large 
condensate fluctuation at ultralow temperatures when 
the ground state is macroscopically populated, the ultra- 
cold ideal gas of atoms has to be described either in the 
microcanonical or the canonical ensemble [55H40j . The 
former one assumes the perfect isolation of the system 
from the environment, while the latter one assumes that 
the system is in contact with a heat bath of certain tem- 
perature: ksT =1/(3. Both ensembles correctly describe 
the fluctuations and correlations of an ideal gas at low 
temperatures. In our approach we apply the canonical 
ensemble. Its partition function Z(j3, N) can be defined 
as 



(£i»i,JV)> 



m— o 7i2— o 



Tlaa— 



where rii denote the number of particles occupying the 
eigenstate with energy e,, and N is the total number 
of particles. The presence of a discrete delta function 
d assures that only partitions with the total number of 
particles equal to N contribute to the sum. For a nonin- 
teracting system the partition function and all the other 
statistical quantities may be computed using the recur- 
rence formulas (see Appendix A for details), based on the 
formula obtained in |41j : 

N 

Z(f3,N) = Y,J2 e ~ nf3e " Z (P> N - n ^ ( 8 ) 

n—1 if 

where we should take Z(j3, 0) = 1. 



C. Ground state population behaviour 
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FIG. 3: Similar to Fig. [2] but for the first excited state. For 
A = 4.5 J the dominating left peak is two orders of magnitude 
stronger than the other. 



Having calculated the partition function, we can get 
the ground state population and its fluctuations. The 
numerical results are presented on Figure [4] First we ob- 
serve the growth of fluctuations. However, as the ground 
state population decreases with temperature, so should 
its fluctuations. The maximum of SNq occurs at cer- 
tain characteristic temperature which depends on disor- 
der strength. We have developed an analytical model to 
explain this behaviour and to give some estimate on the 
characteristic temperature of the maximum of fluctua- 
tions. The shape of the energy spectrum for the lowest 
states may be approximated by a parabola. For an ideal 
gas with parabolic energy spectrum e„ = an 2 all the sta- 
tistical quantities can be calculated analytically. In one 
dimension one can introduce some characteristic temper- 
ature T c = (SaN/it 2 , which determines the regime when 
the ground state becomes macroscopically populated (see 
Appendix B for details). Below T c the ground-state oc- 
cupation number can be calculated using, for instance, 
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FIG. 4: Top: number of ground state atoms obtained in 
canonical ensemble for different values of A, compared with 
the analytical model predictions. Bottom: fluctuations of the 
number of ground state atoms. N = 300. 



the technique of the Maxwell-Demon ensemble 
This yields 



(No) 
N 



1 



T 



(9) 



Above the characteristic temperature too many atoms 
become excited and the Maxwell-Demon method is not 
applicable. We observe that the model is reliable below 
temperature at which about a half of the particles be- 
come excited. The model based on the Maxwell-Demon 
ensemble gives also the correct value of the characteris- 
tic temperature T c at which the fluctuations reach the 
maximum. 



and for various temperatures of the atomic gas. The 
peaks correspond to localized states that are centered at 
various lattice sites. At low temperatures the peaks are 
distributed symmetrically around the central peak, corre- 
sponding to the ground state. The remaining peaks result 
from the contribution of excited states. As the tempera- 
ture increases, the number of populated states gradually 
grows and so does the number of peaks. Similar effect 
can be observed for fluctuations. It turns out that the 
behaviour of fluctuations can be qualitatively understood 
assuming thermal character of the fluctuations for each 
lattice site separately 



(6 2 n i ) = {n i ){{n i ) + l), 



(10) 



Similar result can be also obtained when considering the 
lattice as a set of separated potential wells, and describ- 
ing the statistics of the single well within the grand- 
canonical ensemble in the equilibrium with the rest of 
the lattice sites, which can be treated as a reservoir |43j . 

This approximation 
The com- 



As a result we obtain formula (10). 
works particularly well at high temperatures, 
parison of exact results and the model is presented on 
Figure [6] 



E. Correlations between sites 

For completeness of the analysis we present the corre- 
lations {nifij) — (rii) (rij) between the number of particles 
in different lattice sites, calculated in the canonical en- 
semble. The correlations for various values of disorder 
strength are shown in Fig. [7] It appears that the cor- 
relations between strongly occupied sites are the largest. 
For large disorder the correlations are mainly negative, 
which results from the conservation of the total number 
of particles in the canonical ensemble: the more particles 
occupy certain localized state, the less are left for the 
other states. The positive values appear only between 
the lattice sites where the localization occurs. In con- 
trast, at low disorder, the correlations are positive only 
between neighboring sites, which correspond to the diag- 
onal and to the corners of the graph. This results from 
the fact that at low disorder all the excited states are 
spread along several sites. 



IV. STRONGLY INTERACTING GAS IN A 
QUASI-PERIODIC POTENTIAL 



D. Mean number and fluctuations 

In this section we analyze the mean and fluctuations of 
the number of particles in the wells of the optical lattice. 
In Figs. [5] and [6] we show the mean occupation numbers 
and its fluctuations calculated for some sample parame- 
ters: N — 300 particles, disorder amplitude A = 2.5 J, 



A. Hamiltonian 

When the particles are strongly interacting, we may 
neglect the tunneling term in the Hamiltonian ([!]) in com- 
parison to the remaining terms. This yields 

H = U^2n k (n k - 1) + A^n fc cos(27r/3fc). (11) 

k k 
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FIG. 5: Mean number of particles in each lattice obtained in 
the canonical ensemble for M = 144 sites, N = 300 particles, 
fixed value of the disorder strength A = 2.5 J, calculated for 
various temperatures. 



FIG. 6: Number fluctuations of particles in each lattice site 
obtained in the canonical ensemble for M = 144 sites, N = 
300 particles, fixed value of the disorder strength A = 2.5 J, 
calculated for various temperatures. The insets show com- 
parison of exact results (black) with thermal approximation 
(blue, dashed) around the central peak. 



Now, there are only two energy scales given by U and A, 
and in the subsequent analysis we will express A in units 
of U. In the Hamiltonian (11) the different lattice sites 



are decoupled, thus the only correlation between sites is 
due to the conservation of the total number of particles. 



For such a system the partition function 
Z(N, /3) = V...y e -^ w-W E Mm-D/2 S . 



(12) 

can be calculated exactly using a recurrence relations. 
We have developed a recurrence algorithm to calculate 
the partition function (12), which can be derived by 
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FIG. 7: Correlations between sites {riirij) — (m) {rij), for Aubry-Andre model with ksT = 0.05 J and A = 0.5 J (left); 
k B T = 0.05 J and A = 2.2 J (middle); and for strongly interacting gas with k B T = 0.5(7, M = 89, N = 2M and A = 5U 
(right). 



adding one lattice site in each step of the recurrence (see 
Appendix C for details). 



B. Statistics of the strongly interacting gas in the 
presence of disorder 

We have analyzed the mean particle number and fluc- 
tuations in the wells in the case of strongly interacting 
gas in quasi periodic potentials. As the tunneling process 
is neglected, the localization is not presents and the only 
effects influencing the mean and fluctuations results form 
the variation of the chemical potential at different lattice 
wells. This statement is confirmed by the analysis of 
numerical results shown in Fig. [8] We have performed 
numerical calculations for a moderate-size system con- 
taining M = 89 sites and N = 2M particles. We observe 
that both mean and fluctuations vary stronger from site 
to site as the amplitude of disorder increases. The corre- 
lations between occupation numbers at different sites for 
some example value od A are presented in Fig. [7] As the 
sites are uncoupled in the Hamiltonian (11), the correla- 



tions result only from the constraint on the total number 
of particles, and they are strongest for sites with highest 
occupation numbers. 



V. PROBING THE STATISTICAL PROPERTIES 
OF THE SYSTEM WITH LIGHT SCATTERING 
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FIG. 8: Sample graph for mean occupation number (top) 
and its fluctuations (bottom) for M = 89 sites, N = 2M, 
fesT = 0.5 U and various A. 



We consider the possibility of distinguishing between 
different many-body phases of ultracold atoms in disor- 
dered potentials. This can be done, for instance, by mea- 
suring the properties of the correlation functions. One of 
the possible tools that can bring information about the 
correlation function is the measurement of the proper- 
ties of light scattered on ultracold atoms. Previously, the 
atom-light interactions were suggested as a method to de- 
tect Bose-Einstein condensation in an ultracold gas [25] . 
BCS transition in ultracold fermions [27], statistics of 
ultracold atomic gases or distinguishing between 

quantum phases of ultracold atom in optical lattices |28F 



eh. 

Let us consider the gas of ultracold atoms in an ex- 
ternal potential interacting with a weak and far-detuned 
laser with frequency lu^. Treating the atoms as two-level 
systems, it is possible to adiabatically eliminate the ex- 
cited state and obtain the effective hamiltonian. In this 
way we can calculate the mean number of scattered pho- 
tons with wave vector k and polarization e per unit time 
per solid angle [2"5] 



d 2 N 
dtdn 



n 2 £ ir 



si 



(13) 
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FIG. 9: Experimental setup for light scattering. The optical 
lattice is illuminated by the laser set at angle a. The detector 
is set at angle /3. 



where f2 = Eie^d/h is the Rabi frequency, Sl is the de- 
tuning of the laser, c k = g k ekd/h, q = k k L , |fc| = \k L \ 
(elastic scattering) , Ei stands for the electric field of the 
laser with polarization cl, d is the atomic dipole moment, 
g k is the coupling constant and 



F(q) = d 



(V t (x)V(x)V t (y)V(y)) 



carries the information about the statistics of the system. 
Function F(q) is defined as the Fourier transform of the 
second correlation function. It is equivalent to the static 
structure factor j29j and in the rest of the work we will 
refer to F(q) as the static structure factor. 



In this regime we may also use gaussian approximation 
of the Wannier states. Formula ( 15 ) simplifies to 



El/o(q) 



2 iq(r„-r m ) 



where 



/o(q)= / d 2 r \w (r)\ 2 e 



zqr 



\n n n m ) , (16) 



(17) 



The term q(r n — r m ) may be rewritten as jin — to), 
where 



A 

7 = 7T — (cos p — cos a) , 



(18) 



A is the wavelength of the laser forming the primary opti- 
cal lattice and A/, is the wavelength of the probing laser. 
Angle f3 is the angle at which the detector is set and a is 
the angle of the probing laser (see Figure [9| . 

It is instructive to split F(q) into two parts F c i ass and 
F qU ant [2H] , where 



F c lass(q) = |/o(9) T 



(n m ) 



(19) 



represents the so called classical component of the scat- 
tered light. It is obtained by calculating the average 
|(ikA)| 2 , which is proportional to the amplitude of the 
electric field square. The difference between the total 
function F(q) and the classical part defines the quantum 
component [35] 



A. Light scattering from bosons in an optical 
superlattice 

We now show how to extract the information on cor- 
relations from the intensity of light scattered at different 
angles. For a single atom, one can show that F(q) = 1 
[25] , so F represents the difference between scattering 
from one atom and from the many-body system. In the 
following, we will focus solely on the properties of the 
structure factor. By expanding the field operators into 
Wannier states, we get an equivalent formula for F: 



F(q) = £ <n|e* 

n,n' ,m,m' 



%') (m \e 4qr | to') (gig n ' gLgm') 

(15) 

The matrix elements (n |e lqr | n') are calculated between 
Wannier states localized in sites n and n' . We will con- 
sider the deep lattice regime where Wannier states are 
strongly localized and hence the terms with n ^ n' are 
negligible. This approximation is valid when the lattice 
potential depth is of the order of several recoil energies. 



F quan M) = \ f (q)\ 2 Y, e * q(rn ~ rm) (KO - W (O) 

n , rn 

It gives information about quantum statistical effects in 
the system. Splitting F(q) into these two parts is partic- 
ularly useful when comparing Mott insulator and super- 
fluid phases, as both of them are homogenous so they 
differ only in the quantum component |28l I29| . Here 
this will not be the case, as the system is inhomogeneous 
and already the classical components of various quantum 
phases are different. 

For the homogenous phase with density n, the classical 
part of F can be expressed as 



F c ia SS (q) = |/o(9)| n 



, 2 sin 2 (M7tA/2Al(cos/3 — cos a)) 



sin (7rA/2Ai(cos j3 — cosa)) 



(21) 

which gives us intuition that as M or A/Al increases, 
F(q) should oscillate faster. This quantity has already 
been measured in experiment for a two-dimensional Mott 
insulator |2"H. 




FIG. 10: The averaged spectrum of scattered photons for Aubry-Andre model for different values of A and T. The probing 
beam is set at the angle a = tt/2. For A < 2 J the system is in the superfluid state while for larger A it is in the localized 
regime. 



B. Scattering from localized and delocalized phases 

We now use the method described above to analyze the 
possibility to distinguish localized and delocalized phases 
in Aubry-Andre model. We will use the results obtained 
in the canonical ensemble and presented in the previ- 
ous chapter. There are many parameters which can be 
varied in calculations and in experiment: the number of 
particles N, number of lattice sites M, temperature T, 
primary lattice depth si, probe laser wavelength A and 



the angle at which the detector is set a. We set N = 300, 
M = 144, si = 5, A = A; (the lattice laser wavelength), 
and a = tt/2 and examine how the spectrum of scattered 
photons changes with growing A and temperature. 

As shown on Figure [lOj the growth of A causes ad- 
ditional interference peaks to emerge. This results from 
the influence of the second lattice, which generates addi- 
tional momenta k^, kx — &2 etc. in the system. Similar 
observation was made in [TT], where the impact of the 
second lattice on the noise correlations was studied ex- 
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perimentally. As A crosses the transition point, due to 
incommensurability of the lattices, the angular distribu- 
tion flattens as the number of interference peaks goes to 
infinity. As a result we are able to detect localization 
for high A, as well as observe the growing impact of the 
secondary lattice for low disorder. 

High temperature rises the number of excited particles 
and disturbs the angular distribution of photons in two 
ways. Below the transition point A = 2 J it reduces the 
visibility of the interference peaks. For higher A the pres- 
ence of several localized states produces the interference 
peaks in the distribution which would not be present at 
T = 0. 



A = 




C. Scattering from strongly interacting gas 

We now analyze the scattered spectrum for the 
strongly interacting gas, keeping the parameters of the 
bichromatic lattice unchanged and setting M = 89 and 
N = 2M. Sample pictures are shown on Figure [TT] In 
the absence of disorder A = 0, the ultracold Bose gas 
forms a Mott insulator, and in such a case it was pre- 
dicted that the angular distribution of scattered photons 
should exhibit the pattern of interference fringes, with a 
set of minima where there are no scattered photons [35] . 
This is due to the contribution from the classical compo- 
nent, while the quantum part is zero due to the absence 
of the correlations. In contrast to the previous works, in 
our calculations we include the effects of the correlations 
between wells due to the constraint on the total number 
of particles in the canonical ensemble. This gives rise to 
the nonzero quantum component, and as the result the 
minima are no longer present in the angular distribution. 

For growing disorder we again observe the appearance 
of new peaks in the spectrum. In this case they reflect 
the fact that the additional lattice of different period was 
added and the gas has a new density profile. However, 
there is no localization, so there is no qualitative change 
for growing A. The angular distribution will flatten only 
in the limit A — > oo. This means that in principle it 
should be possible to distinguish the strongly interacting 
phase from the localized phase which may be useful in 
examining the role of interactions in disordered systems. 
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FIG. 11: The angular distribution of scattered photons for 
the strongly interacting gas for different values of A and T. 
The probing beam is set at the angle a = 7r/2. 



VI. CONCLUSIONS 

In conclusion, we studied the statistical properties of 
a Bose gas confined in a bichromatic optical lattice at 
finite temperature. We considered two limits, when the 
Hamiltonian can be diagonalized exactly: the ideal gas, 
when there is no interactions, and strongly interacting 
gas, when one can neglect the inter-well tunneling. We 
analized the mean, fluctuations and correlations between 
lattice sites occupation numbers for the Bose-condensed 
phase, localized phase and strongly interacting phase. 



We have shown that some important information about 
the structure factor can be extracted using light scatter- 
ing, which makes possible to distinguish different phases 
and explore the phase diagram experimentally. 
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Appendix A: Calculating the statistical quantities 
for an ideal gas 

In this appendix we derive a recurrence formula allow- 
ing for calculation of the partition function of the nonin- 
teracting gas exactly [H]. Let us denote the probability 
of finding exactly n particles in a given state p as p^ (n) . 
It may be calculated as p^ (n) —p^ (n+1) which represent 
the probabilities of finding at least n (or n+1, respec- 
tively) particles in an eigenstate p. Probability p^ (n) is 



given by the formula 



P " W Z(J3,N) ^ 



E ■•■ E e-^^^i 



,N)- 

(Al) 

By changing the summation index fi^ = — n we ob- 
tain exactly a formula for Z(J3, N — n) with prefactor 
e - n P e »/Z(P,N). Therefore, 



(n) = e -ngc„ m N-n) c _ (n+1)fUu m N-n-l) 



Z((3,N) 



(A2) 

With this result we can easily calculate the mean oc- 
cupation number of state p, directly from its definition 



= J2n=i n Pfj.( n ) an( i using ( A2), obtaining 



1 N 

w = y- E 



e ^ N—n • 



(A3) 



71=1 



Summing (n^) over all states, we obtain the desired for- 
mula for the partition function 



N 



Z(/3,A0 = EE' 



,-nPe u 



Z(fi,N-n). 



(A4) 



n=l 



For calculating the fluctuations of the occupation number 
for a single site, as well as the correlations of the occupa- 
tion numbers between different sites, we need to express 
(nkni) in terms of the partition function. Let us denote 
the probability of finding exactly n particles in state a 
and m particles in state 7 by P a-,(n , to). This parameter 
may be calculated similarly to (A2), using p^ (n,m) de- 
fined as the probability of finding at least n particles in 
state a and at least m in state 7. This quantity fulfills 
the following relation: 

p a7 (n, m) =p=^ 7 (n, m) + p^ 7 (n + 1, m + 1) - 

-p| 7 (n, m + 1) -p^ 1 {n + l,m) (A5) 

Similar calculations as for the mean occupation yields 

p%{n, to) = e -«*«-«*e, m z {p ™~ U) (A6) 



Next, we calculate (n Q n 7 ) from the definition: (n a n 7 ) — 
J2n m nm Pai( n ' m ) an d usm g (A5|. After straightfor- 
ward calculations, we get 



1 



N N 

Z(B,N) ^ ^ ' 

' n =l m=l 



-n/3e a — m/3e- 



(A7) 

Now, we are in position to calculate the mean occupa- 
tion numbers and correlations between different sites of 
the optical lattices. Formally, correlations between two 
lattice sites are given by the trace of density matrix p 
with operators hifij. For the canonical ensemble p = 

e-P"/Z(P,N), so {riiUj) = Z-^^N^r^e-^hinjY 

Because p has a simple form only in the basis of hamilto- 
nian eigenstates, we have to find the relation between cre- 
ation and annihilation operators of single-particle eigen- 
states b a and lattice sites §i numerically. By expressing 
the mean values of site operators by mean values of eigen- 
state operators and utilizing the fact that the number of 
atoms is constant so only certain terms are non-zero, we 
obtain 



: K) + Ei c ° 



\c]\ (n Q n^)- 



+ <t + 2^ c r^cfc^ (n a ) 

= E l c "! 2 \ c l\ 2 K«i>+E cr4cT^{{n a n v ) + (n Q )), 

(A8) 

as well as the formula for the mean occupation numbers 



Ei c "i 2 <^> 



(A9) 



Appendix B: Ideal Bose gas with quadratic energy 
spectrum 

Let us consider a gas of noninteracting bosons in an 
external trap with the quadratic energy spectrum: efe = 
ak 2 , where a is constant. The logarithm of the grand 
canonical partition function is by definition given by 



00 

ln3(z,)8)= J>(l 



-0ak' 



(Bl) 



fe=0 



where z = e ^ M is called the fugacity. The expectation 

value of the total number of particles is given by (N) — 

d 

Oz 



z-g- kS(z, (3), which yields 



(N) = 



^ ze~^ k 



E 



P(p,-ak 2 )/2 



Z^l- Z e-^ 2 ^ 2sinh((/3afc 2 - ftp) /2)' 

ft — k — 

(B2) 
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In the low-temperature limit: <C 1 and (3 a <C 1, we 
may keep the first term in series expansion of sinh() and 
exp functions, obtaining 



(N) 



1 - z 



+ 



1 



a/3 ^ k 2 
r fe=i 



1 - z 



^ <"»> 



The ground state occupation (No) is equal to z/(l — z). 
At the characteristic temperature T c it becomes macro- 
scopic. From this observation we conclude that the 
characteristic temperature of the gas is given by T c = 
6aN/-K 2 , and that below T c 



(No) 
(N) 



1 



T 

r. 



(B4) 



tonian separates the lattice sites and we can develop a 
different method. We start from definition 



ni riM 

(CI) 

Due to the separation of the lattice sites the partition 
function for the first m lattice sites Z(N, (3, m) can be 
expressed by Z(N, /3,m— 1): 



Appendix C: The partition function for the strongly 
interacting gas 

For the strongly interacting gas the previous method 
does not work but we may use the fact that the hamil- 



Z N (p, m) = Y J e-^ n --^ n --^Z N - nm (ft m - 1). 

(C2) 

Calculation of the occupation numbers, fluctuations and 
correlations may be done in a similar manner. 



[1] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[2] H. Hefei, A. Strybulevych, J. H. Page, S. E. Skipetrov, 

and B. A. van Tiggelen, Nat. Phys. 4 (2008). 
[3] A. A. Chabanov, M. Stoytchev, and A. Z. Genack, Nature 

404, 850 (2000). 
[4] D. S. Wiersma, P. Bartolini, A. Lagendijk, and R. Righ- 

ini, Nature 390, 671 (1997). 
[5] G. Roati, C. D'Errico, L. Fallani, M. Fattori, C. Fort, 

M. Zaccanti, G. Modugno, M. Modugno, and M. Ingus- 

cio, Nature 453, 895 (2008). 
[6] J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, 

P. Lugan, D. Clement, L. Sanchez-Palencia, P. Bouyer, 

and A. Aspect, Nature 453, 891 (2008). 
[7] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
[8] D. Clement, A. F. Varon, M. Hugbart, J. A. Retter, 

P. Bouyer, L. Sanchez-Palencia, D. M. Gangardt, G. V. 

Shlyapnikov, and A. Aspect, Phys. Rev. Lett. 95, 170409 

(2005). 

[9] D. Clement, A. F. Varen, J. A. Retter, L. Sanchez- 
Palencia, A. Aspect, and P. Bouyer, New Journal of 
Physics 8, 165 (2006). 
[10] J. E. Lye, L. Fallani, C. Fort, V. Guarrera, M. Modugno, 
D. S. Wiersma, and M. Inguscio, Phys. Rev. A 75, 061603 
(2007). 

[11] V. Guarrera, N. Fabbri, L. Fallani, C. Fort, K. M. R. 
van der Stam, and M. Inguscio, Phys. Rev. Lett. 100, 
250403 (2008). 

[12] D. Thouless, Physics Reports 13, 93 (1974). 

[13] P. W. Anderson, D. J. Thouless, E. Abrahams, and D. S. 
Fisher, Phys. Rev. B 22, 3519 (1980). 

[14] H. von Dreifus and A. Klein, Communications in Math- 
ematical Physics 124, 285 (1989), 10.1007/BF01219198. 

[15] S. Aubry and G. Andre, Ann. Isr. Phys. Soc. 3, 33 (1980). 

[16] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 



Fisher, Phys. Rev. B 40, 546 (1989). 
[17] B. Damski, J. Zakrzewski, L. Santos, P. Zoller, and 

M. Lewenstein, Phys. Rev. Lett. 91, 080403 (2003). 
[18] G. Roux, T. Barthel, I. P. McCulloch, C. Kollath, 

U. Schollwock, and T. Giamarchi, Phys. Rev. A 78, 

023628 (2008). 

[19] L. Fontanesi, M. Wouters, and V. Savona, Phys. Rev. A 

81, 053603 (2010). 

[20] M. Greiner, O. Mandel, T. Esslinger, T. Hansch, and 

I. Bloch, Nature 415, 39 (2002). 
[21] R. Roth and K. Burnett, Phys. Rev. A 68, 023604 (2003). 
[22] S. Foiling, F. Gerbier, A. Widera, O. Mandel, T. Gericke, 

and I. Bloch, Nature 453, 481 (2005). 
[23] E. Altman, E. Demler, and M. D. Lukin, Phys. Rev. A 

70, 013603 (2004). 
[24] C. Weitenberg, P. Schaufi, T. Fukuhara, M. Cheneau, 

M. Endres, I. Bloch, and S. Kuhr, Phys. Rev. Lett. 106, 

215301 (2011). 

[25] M. Lewenstein and L. You, Phys. Rev. Lett. 71, 1339 
(1993). 

[26] Z. Idziaszek, K. Rzazewski, and M. Lewenstein, Phys. 

Rev. A 61, 053608 (2000). 
[27] W. Zhang, C. A. Sackett, and R. G. Hulet, Phys. Rev. A 

60, 504 (1999). 

[28] I. Mekhov, C. Maschler, and H. Ritsch, Nat. Phys. 3, 319 
(2007). 

[29] K. Lakomy, Z. Idziaszek, and M. Trippenbach, Phys. Rev. 

A 80, 043404 (2009). 
[30] S. Rist, C. Menotti, and G. Morigi, Phys. Rev. A 81, 

013404 (2010). 

[31] J. S. Douglas and K. Burnett, Phys. Rev. A 84, 033637 
(2011). 

[32] J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper- 
Kurn, D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 

82, 4569 (1999). 



12 



[33] K. Eckert, O. Romero-Isart, M. Rodriguez, M. Lewen- 

stein, E. Polzik, and A. Sanpera, Nat. Phys. 4, 50 (2008). 
[34] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[35] M. Modugno, New Journal of Physics 11, 033023 (2009). 
[36] C. Aulbach, A. Wobst, G.-L. Ingold, P. Hnggi, and 

I. Varga, New Journal of Physics 6, 70 (2004). 
[37] S. Y. Jitomirskaya, The Annals of Mathematics 150, pp. 

1159 (1999). 

[38] H. D. Politzer, Phys. Rev. A 54, 5048 (1996). 



[39] M. Gajda and K. Rzazewski, Phys. Rev. Lett. 78, 2686 
(1997). 

[40] P. Navez, D. Bitouk, M. Gajda, Z. Idziaszek, and 
K. Rzazewski, Phys. Rev. Lett. 79, 1789 (1997). 

[41] C. Weiss and M. Wilkens, Opt. Express 1, 272 (1997). 

[42] S. Grossmann and M. Holthaus, Opt. Express 1, 262 
(1997). 

[43] Z. Idziaszek, Ph.D. thesis, Center for Theoretical Physics, 
Polish Academy of Sciences (2002). 



